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摘要 : 光 干 涉 技术 需要 对 光 程 差 进行 精确 的 补偿 和 控制 ， 因 此 光 程 差 的 精确 探测 尤为 重 
要 。 主 要 针对 色散 条 纹 法 在 光 程 差 探 测 中 的 应 用 进行 讨论 ， 通 过 理论 分 析 、 仿 真 计 算 和 实 
验 ， 验 证 了 在 一 个 波长 范围 内 的 光 程 差 与 色散 干涉 条 纹 纵向 光 强 峰值 偏 移 量 之 间 存 在 线性 关 
系 。 分 析 实 验 获得 的 白光 色散 干涉 条 纹 数据 发 现 ， 理 论 计 算 预测 值 有 一 定 的 偏差 ， 通 过 对 实 
验 系统 中 各 主要 误差 源 进行 分 析 发 现 ， 分 光 棱 镜 在 运动 过 程 中 引入 的 俯仰 起 伏 是 导致 实验 误 
差 的 主要 原因 。 
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光学 综合 孔径 技术 ， 又 称 光 学 合成 孔径 技术 ， 是 将 多 个 较 小 口径 的 光学 元 件 或 光学 系统 进行 精确 
的 排列 ， 使 通过 各 子 孔径 的 光束 在 焦 平 面 上 满足 一 定 的 同 相 位 条 件 ， 实 现 相干 释 加 ， 经 图 像 复 原 处 理 
后 ， 达 到 与 单一 较 大 通 光 口径 相当 的 衍射 极限 分 辨 率 的 技术 。 

H. Fizeau 最 早 在 1868 年 提出 恒星 光 干 涉 技术 ,希望 利用 干涉 法 测量 恒星 的 角 直 径 。1874 年 ， 
Steven 通过 在 Large Foucault 望远镜 前 增加 掩 膜 的 方法 ， 首 次 成 功 观 测 到 恒星 的 干涉 条 纹 。 到 目前 为 
ik, 世界 上 已 经 建立 了 多 个 拼接 子 望 远 镜 的 光学 综合 孔径 实验 系统 。 美 国 麻 省 理工 学 院 ( Massachusetts 
Institute of Technology ，MIT)2002 年 研制 的 GOLAY-3 型 光学 卫星 地 面 样机 (ARGOS)'" ， 演 示 验 证 了 空 
间 光 学 综合 孔径 技术 的 可 行 性 ， 其 等 效 孔径 为 0.62 m， 角 分 辨 率 为 0.35”( 优 于 单个 子 孔径 的 角 分 辨 
率 0.55") ， 视 场 角 约 为 3'。 美 国 洛克 -马丁 公司 研制 的 Staro 装置 是 由 9 个 子 孔径 组 成 的 Fizeau 型 
光学 综合 孔径 望远镜 ， 单 子 孔径 为 125 mm， 经 过 共 相 位 干涉 ， 实 现 了 与 610 mm 口径 衍射 极限 相当 的 
分 辨 率 ， 视 场 大 小 为 1 urad。 

北京 理工 大 学 于 2009 年 成 功 研 制 了 共 相 位 检测 与 控制 平台 ,对 拼接 镜面 成 像 系统 闭环 控制 ,最 
终 形成 干涉 条 纹 '。 我 国 自主 研制 的 LAMOST 望远镜 于 2009 年 竣工 验收 ， 采 用 主动 光学 技术 ， 有 34 
块 拼接 主 镜 和 28 块 拼接 施 密 特 校 正 板 ,但 是 只 实现 了 共 焦 成 像 ， 分 辩 率 仍然 受 限 于 单子 镜 的 口径 。 中 
国 科 学 院 西 安 光 学 精密 机 械 研 究 所 于 2011 年 进行 了 综合 孔径 实验 系统 的 搭建 ， 开 展 了 理论 分 析 、 仿 真 
研究 和 有 限 距离 综合 孔径 原理 验证 实验 。 上 海天 文 台 也 于 2011 年 建成 2 个 望远镜 组 成 的 综合 孔径 成 像 
系统 ， 采 用 棱锥 波 前 探测 ， 但 未 实现 共 相 成 像 。 与 国外 相 比 ， 我 国 的 天 文 光 干涉 技术 有 竺 进一步 发 展 。 

光学 综合 孔径 望远镜 各 子 孔径 间 达 到 共 相 时 才能 形成 干涉 ， 实 现 复 振幅 的 全 加 ， 提 高 空间 分 辨 
率 ， 和 否则 ， 只 是 光 强 至 加 ， 系 统 分 辩 率 仍 受 限 于 子 孔 径 的 衍射 极限 。 因 此 ， 各 光束 间 光 程 差 的 精确 探 
测 对 光学 综合 孔径 望远镜 光 干 涉 的 实现 至 关 重 要 。 

目前 国内 外 发 展 了 多 种 光 程 探测 的 方法 。 离 焦 面 检测 : (1) 曲率 传 感 法 站 通过 测量 子 镜 边 缘 的 相 
位 不 连续 性 重 构 子 镜 间 的 光 程 差 误 差 ，(2) 相 位 恢复 法 中 通过 焦 面 和 一 对 离 焦 面 的 光 强 数据 反 解 瞳 面 
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相位 分 布 计算 光 程 差 误 差 。 出 瞳 面 检测 : 四 棱锥 波 前 检测 法 ”通过 检测 4 个 出 瞳 像 点 能 量 差异 计算 光 
程 差 误差 。 焦 面 检测 法 : (1) 加 州 大 学 Gary Chanan 提出 的 宽带 哈 特 曼 - 夏 克 法 和 窄带 哈 特 曼 - 夏 克 法 ， 
通过 采样 孔径 的 衍射 图 样 与 样板 图 样 的 相关 系数 计算 拼接 子 镜 间 的 光 程 差 误差 ; (2) 和 詹姆斯 韦伯 太 
空 望远镜 使 用 了 色散 条 纹 法 '“"， 通 过 处 理 相 邻 子 孔径 间 色 散 干 涉 条 纹 ， 对 光 强 分 布 进行 拟 合 ， 得 到 
光 程 差 误差 。 

国家 天 文 台 承 担 了 863 计划 光学 稀 玻 综合 成 像 系 统 原 理 样 机 和 综合 数据 处 理 技术 的 课题 ， 本 文 以 
此 为 背景 ， 对 光 程 差 探测 方法 开展 了 理论 和 实验 研究 。 针 对 系统 的 情况 ， 选 择 色散 条 纹 法 探测 一 个 波 
长 范围 内 的 光 程 差 。 在 分 析 双 孔 干 涉 理论 和 仿真 模拟 的 基础 上 ， 讨 论 了 光 程 差 与 干涉 条 纹 光 强 纵向 峰 
值 偏 移 量 之 间 的 关系 ， 搭 建 了 实验 装置 ， 通 过 处 理 干涉 条 纹 得 到 了 光 程 差 ， 并 对 实验 误差 进行 了 系统 
分 析 ， 找 到 了 影响 计算 精度 的 主要 原因 是 运动 平台 的 俯仰 起 伏 ， 为 后 续 改 进 指明 了 方向 。 


1 理论 基础 


双 圆 孔 在 望远镜 焦 平面 处 ， 经 沿 x 方向 色散 的 AMICI 棱镜 ， 形 成 的 干涉 图 样 强 度 分 布 ”为 


(ma? [25,(Z) 7? 2m (dy J| 
re 9 =2(55) | Z | 个 区 上 


其 中 ，2Z=2mwar/Af; a 为 圆 孔 半 径 ; r= Vx ty’; 1/ 为 望远镜 焦距 ; d 为 基线 长 度 ; Ia, V 为 色散 条 纹 
信和 号 光 强 值 ; /为 色散 条 纹 信号 振幅 ; A (x) 为 沿 x 方 向 色散 的 波长 ; 4$ 为 两 圆 孔 干 涉 的 光 程 差 ; x，y 


二。 分 别 为 横 纵 坐标 。 

Fizeau 型 光 干 涉 系 统 中 光 程 差 探 测 部 分 双 圆 和 孔 示 意图 如 图 1， 其 中 参数 为 : 圆 孔 半径 a=10 mm; 
< 一 基线 长 度 d=40 mm; 焦距 f= 1 680 mm, 

~ d= 40.00 mm 

N a= 10.00 mm 


图 1 两 圆 孔 参 数 


Fig. 1 The parameter of two circle holes 


先 讨论 单 波长 下 的 光 干 涉 ， 令 x=0， 则 Z=2moyAf， 此 时 只 需要 对 干涉 条 纹 纵向 (y 方向 ) 的 光 


强 分 布 进行 分 析 : 
2J (Pret) 
(ma)? VF 2m (dy 
50 =2(55) 2may [i seos [5 (Fra) |} (5 
Mf 
在 -AX2~ +AZ2 之 间 ， 自 变量 光 程 差 $4 与 因 变 量 峰 值 偏 移 量 Yaa LIA A CHP a WEEER ) : 
Af Fina A 
zdy Brees ( 2 =) a~ 1) = = 7 Q. (3) 
2 J, 
Af 
2TA Y rox 


Af 
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2 色散 条 纹 法 仿真 模拟 


假设 自 变量 光 程 差 $ 在 一 个 波长 范围 内 ， 此 时 由 与 峰值 偏 移 量 yw 之 间 存 在 线性 关系 ; 
E K (常数 ) (4) 


将 系统 参数 代入 (4) 式 ， 分别 计算 500 nm( ZL) , 600 nm( 绿 )、700 nm( 蓝 ) 3 种 单 波长 干涉 条 纹 
纵向 光 强 峰值 偏 移 量 Yal PEH y 代替 ) 与 光 程 差 (PIAA piston 代替 ) 之 间 的 关系 。 利 用 MATLAB 
进行 数值 模拟 ， 计 算 结果 如 图 2， 其 中 横 坐 标 光 程 差 的 采样 间隔 为 10 nm， 纵 坐标 y 为 各 光 程 差 下 干 
涉 条 纹 纵向 光 强 对 应 的 峰值 位 置 。 取 线性 相关 系数 rE CUT : 


By (5) 


E 1 n E 1 n n = n 
其 中 , x= 一 Xis 7= 一 > Vis L=}, (x; =x)’; l, = 
n n i=l 


之 之 Le HN G by = a 7)" 

在 3 种 单 波长 条 件 下 ,计算 K 和 7 值 ， 均 得 到 KK 二 -39. 51， 线 性 相关 系数 > 优 于 99. 99%。 因 此 ， 可 
认为 在 单 色光 干涉 时 ， 光 程 差 在 一 个 波长 范围 内 与 y 满足 线性 关系 ， 并 且 不 同 波长 对 斜率 的 影响 较 小 。 

对 多 波长 色散 干涉 条 纹 进行 分 析 ， 假 设 线 色散 关系 x=2 x 10°(A-650 x 10”) (单位 : m), WK 
A=600~700 nm, 采样 间 隔 AA=1 nm， 光 程 差 piston=-250 nm~250 nm 范围 内 ， 得 到 光 程 差 与 峰值 偏 
移 量 y 的 拟 合 曲线 如 图 3， 直线 拟 合 得 到 斜率 KK 二 -39. 37(K=y/piston) ， 线 性 相关 系数 + 二 99. 999 9% 。 
因此 多 波长 光 形 成 的 色散 干涉 条 纹 纵向 光 强 峰值 位 置 y 与 光 程 差 也 呈 线 性 关系 ， 和 斜率 与 单 波长 情况 基 
本 相同 。 
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图 2 峰值 偏 移 量 与 光 程 差 之 间 的 关系 (3 种 波长 ) 图 3 峰值 偏 移 量 与 光 程 差 之 间 的 关系 (多 波长 ) 
Fig.2 Relationship between y and piston Fig.3 Relationship between y and piston 


( single-wavelength ) ( multi-wavelength ) 


3 实验 和 数据 处 理 


在 实验 室 搭建 实验 装置 如 图 4。 实验 系统 由 5 部 分 组 成 : (1) NKT 白光 光源 ; (2) 光 源 望 远 镜 和 
合 光 成 像 望 远 镜 ; (3) CREE S,, S,; (4) AMICI 棱镜 ; (5)CMOS 相机 (2 048 x 1088 像 元 ) 。 系 统 
光路 示意 图 见 图 5。 分 光 棱镜 5, 位 置 保持 不 变 ，5, 安 装 在 直线 运动 导轨 (LPS-45) 上 ， 可 通过 计算 机 控 
制 8 运动 改变 光 程 差 。 实 验 采集 的 零 光 程 干涉 条 纹 图 像 如 图 6。 
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图 4 实验 室 系 统 搭建 照片 


Fig.4 Picture of experimental system 


根据 y 与 光 程 差 线性 相关 的 仿真 结果 ， 在 
-250 nm~ +250 nm 范围 内 ， 以 10 nm 为 采样 间 
隔 ， 分 别 采 集 不 同 光 程 差 下 的 色散 干涉 条 纹 ， 并 
进行 数据 处 理 。 为 了 减 小 误差 ， 对 各 干涉 条 纹 均 
取 n 组 纵向 (垂直 于 色散 方向 ) 光 强 数 据 ， 分 别 
进行 拟 合 得 到 峰值 位 置 偏 移 量 y,，y,,，…，y,， 
求 其 均值 得 到 yo 

图 7 为 实验 中 一 组 干涉 条 纹 纵向 光 强 峰值 偏 
移 量 y 与 光 程 差 的 拟 合 结果 。 其 中 ， 横 坐标 光 程 
差 的 范围 是 -250 nm~ +250 nm， 纵 坐标 y 为 干涉 
条 纹 纵向 光 强 峰值 偏 移 量 的 均值 。 线 性 拟 合 得 到 
的 线性 相关 系数 r=99. 81%, FER K=-50. 74。 


| | 上方 分 光 棱镜 S | 
一 上 下方 分 光 棱 倘 S。 | 
S> 运 动 方向 


CMOS 相 机 CJ 


图 5 系统 光路 示意 图 
Fig.5 Sketch map of system 


图 6 零 光 程 干涉 条 纹 


Fig.6 Interferometric fringe of piston zero 


通过 上 述 结果 可 以 看 出 ， 实 验 数据 处 理 与 仿 
真 计算 得 到 的 斜率 之 间 存 在 一 定 程度 的 偏差 。 接 
下 来 对 系统 的 主要 误差 源 进 行 分 析 。 


4 误差 分 析 


4.1 系统 主要 参数 对 斜率 的 影响 

为 了 分 析 系 统 各 主要 设计 参数 对 拟 合 斜率 
的 影响 ， 利 用 MATLAB 建立 多 波长 色散 干涉 模型 ， 
分 别 分 析 了 合 光 成 像 望 远 镜 焦距 人 圆 孔 半径 a 
和 基线 长 度 d 的 变化 对 斜率 天 的 影响 。 光 源 波长 
范围 设 定 在 600~700 am， 结果 见 表 1、 表 2。 


x10% 


-2 -1 0 1 2 
Piston /m x107 
图 7 干涉 条 纹 峰 值 偏 移 量 与 光 程 差 拟 合 结果 (实验 结果 ) 
Fig.7 Relationship between y and piston 


(result of the experiment) 
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从 表 1、 表 2 可 以 看 出 ， 合 光 成 像 望 远 镜 焦距 /在 + 100 mm 误差 范围 内 ， 造 成 斜率 变化 范围 在 + 
2 以 内 ; ALAA a 在 +1 mm 误差 范围 内 ， 造 成 斜率 变化 范围 约 + 上 0.5 以 内 ; 基线 长 度 d 在 + 上 1mnm ik 
差 范 围 内 ， 造 成 斜率 变化 范围 二 1 以内。 可见， 以 上 3 个 系统 参数 对 斜率 的 影响 较 小 ， 不 是 实验 中 影 
响 斜 率 的 主要 因素 。 


R1 aa 为 10 mm 时, f 和 4d 对 kK 的 影响 表 2 f 为 1680 mm f}, a Md 对 KK 的 影响 
Table 1 The influence of f and d to K, Table 2 The influence of a and d to K, 
When a= 10mm When f= 1680mm 
d 39 mm 40 mm 41 mm d 39 mm 40 mm 41 mm 
f=1580 mm —37. 86 -37. 04 -36.25 a=9 mm -40.76 —39. 86 -39. 00 
f=1 680 mm —40. 24 -39. 37 —38. 54 a=10 mm —40. 24 -39. 37 -38. 54 
f=1780 mm -42. 62 -41. 70 —40. 82 a=11 mm -39. 69 —38. 86 -38. 05 


4.2 光源 光谱 曲线 对 斜率 的 影响 

NKT 白光 光源 的 光谱 曲线 中 ， 各 个 波长 的 光 强 度 值 均 有 所 差异 。 为 了 模拟 光谱 曲线 对 色散 条 纹 
的 影响 ， 假 设 了 两 组 光谱 曲线 : (1) 各 波长 的 光 强 满足 均匀 随机 分 配 ( 以 0.5~1 之 间 的 随机 数 为 比例 
进行 分 配 ) ，(2) 各 波长 的 光 强 按 一 定 规律 进行 分 配 (光谱 曲线 见 图 8(a) ， 仿 真 计算 结果 如 图 8(b) )。 
仿真 结果 显示 : 两 组 光谱 曲线 下 得 到 的 y 与 光 程 差 之 间 的 斜率 玉 均 为 -39.37， 线 性 相关 系数 > 为 
99. 999 9% ， 并 且 与 各 波长 光 强 相同 时 表现 一 致 。 因 此 ， 光 谱 曲 线 也 不 是 造成 料 率 变 化 的 原因 。 
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(a) (b) 
图 8 (a) 光 谱 曲 线 ; (b) 仿 真 结果 的 影响 


Fig.8 (b) The simulation results under the influence of (a) spectral curve 


4.3 ”分 光 棱 镜 运 动 过 程 中 引入 的 角度 扰动 对 斜率 的 影响 

带动 分 光 棱 镜 $ 运动 的 压 电 陶瓷 线性 运动 平台 是 PI 公司 的 LPS-45， 其 最 大 行程 13 mm(+ 6.5 
mm ) ， 双 向 可 重复 性 为 + 18 nm， 最 小 运动 增 量 为 6nm， 角 度 扰 动 (pitch&yaw ) 为 + 50 urad。 在 实验 系 
统 中 ,俯仰 角度 扰动 使 条 纹 发 生 上 下 偏 移 ， 偏 航 角度 扰动 使 条 纹 发 生 左 右 偏 移 。 首 先 ， 对 两 束 光 之 间 
存在 左右 、 上 下 错位 的 情况 进行 仿真 ， 结 有 果 显 示 : 

(1) 两 束 光 左右 错位 (Xshifh) -4 pm~ +4 pm: 光 程 差 为 0 时， 色散 干涉 条 纹 纵向 光 强 峰值 偏 移 量 
的 最 大 值 约 为 0.01 nm; 光 程 差 为 + 100 nm 时 ， 与 未 发 生 左右 错位 情况 下 色散 干涉 条 纹 纵向 光 强 峰值 
偏 移 量 之 差 的 最 大 值 约 为 0.1 nm。 可 见 ， 左 右 错 位 对 色散 干涉 条 纹 光 强 峰 值 偏 移 量 的 影响 较 小 。 
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(2) 两 束 光 上 下 错位 (Ysh 访 )-4 pm~ + 4 pm: 光 程 差 为 0 时， 色散 干涉 条 纹 纵向 光 强 峰值 位 置 偏 


移 量 随 着 Yshif Wt 


兽 加 而 二 


首 加 ， 偏 移 量 最 大 值 约 为 122 nm; 光 程 差 为 + 100 nm 时 ， 色 散 干 涉 条 纹 纵向 


光 强 峰值 位 置 偏 移 量 随 着 Yshifi 的 增加 而 增加 ， 与 未 发 生 上 下 错位 情况 下 色散 干涉 条 纹 纵向 光 强 峰值 


偏 移 量 之 差 的 最 大 值 约 为 126 nm。 可 见 ， 上 下 错位 对 色散 干涉 条 纹 光 强 峰值 偏 移 量 的 影响 较 大 。 
假设 两 色散 条 纹 在 0 光 程 差 下 相互 重合 ， 通 


过 5, 运动 产生 + 10 nm 光 程 差 ， 同 时 引入 上 下 


表 3 上 下 条 纹 错 位 对 斜率 的 影响 
Table 3 K under the influence of two 


错位 -4 hm~ + 4pm, RMR KE KERK, 
为 -27. 20~ -51. 45( 如 表 3) 。 实 验 结果 表明 ， 在 


fringes up and down 


运动 平台 最 大 行程 内 移动 5, 将 导致 探测 器 上 的 上 下 错位 /pm ARK 
ERRARE ETNI l= + 27.75 像素 ， 对 应 (上 移 为 正 ,下 移 为 负 ) 
光束 角度 变化 a=15.5 pm/f= + 90.8 urad, S, f -4 -27. 20 
度 变 化 ao/2= + 45. 4 urad (Æ + 50 urad A), SK -2 -33.25 
统 标 称 的 数据 相符 。 因 此 ， 可 以 判断 上 下 错位 是 P 255 49 
造成 实验 数据 与 仿真 结果 斜率 值 之 间 存 在 较 大 > Per 
偏差 的 主要 原因 ， 这 是 由 运动 平台 的 俯仰 角度 误 age 
差 决 定 的 。 
5 总 结 

文中 主要 介绍 了 使 用 色散 条 纹 法 对 光 程 差 进行 探测 的 方法 和 实验 结果 。 首 先 ， 数 值 计 算 显 示 ， 


光 程 差 在 一 个 波长 范围 内 与 色散 干涉 条 纹 纵 向 光 强 峰值 偏 移 量 之 间 存 在 良好 的 线性 关系 。 其 次 ， 
搭建 实验 光路 ， 开 展 实验 研究 ， 通 过 控制 分 光 棱 镜 5, 运 动 改变 两 光束 间 的 光 程 差 .再 对 采集 的 色散 
干涉 条 纹 进行 数据 处 理 得 到 斜率 数据 ， 发 现 与 仿真 得 到 的 斜率 值 有 较 大 偏差 。 经 过 较 全 面 的 误差 分 
析 ， 分 别 对 主要 系统 参数 、 光 源 光 谱 曲 线 和 分 光 棱 镜 运 动 引 入 的 角度 扰动 进行 了 分 析 ， 结 果 显示 分 光 
棱镜 运动 过 程 中 引入 的 俯仰 角度 误差 导致 的 条 纹 上 下 偏 移 是 实验 数据 与 仿真 计算 斜率 值 偏差 较 大 的 主 


针对 这 一 问题 ,希望 用 系统 标定 的 方式 进行 修正 。 即 利用 双 频 激光 位 移 干涉 仪 对 运动 棱镜 位 置 进 
行 精确 测量 (分 辩 率 达 1 nm) ， 以 此 对 杰 镜 位 置 进行 标定 ， 进 而 得 到 经 标定 的 糙 率 K,。 在 实验 中 ， 用 
标定 的 K, 对 光 程 差 进行 求解 : piston=y/Ko。 
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Abstract; For optical interference, piston error detection influences the accuracy of piston compensation 
considerably. We build an optical experimental system and choose dispersed fringe sensing (DFS) as piston 
detection method. Through numerical simulations and experiments, we find the piston within a wavelength has 
a linear relationship with longitudinal intensity peak offset of dispersed interference fringes. In our experiment , 
we use white light source (100nm spectral width) to generate dispersed interference fringes. The piston will 
be introduced intentionally by moving beam splitter platform. We find there is a difference between simulation 
results and experimental results. Then, we analyze the influence of some experimental parameters on piston 
detection accuracy. Finally, we find the pitch movement disturbance of beam splitter platform is the main 
reason of the difference. 
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